Numerical relativity for D dimensional space-times: 
Head-on collisions of black holes and gravitational wave extraction 
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Higher dimensional black holes play an exciting role in fundamental physics, such as high energy 
physics. In this paper, we use the formalism and numerical code reported in [J] to study the head-on 
collision of two black holes. For this purpose we provide a detailed treatment of gravitational wave 
extraction in generic D dimensional space-times, which uses the Kodama-Ishibashi formalism. For 
the first time, we present the results of numerical simulations of the head-on collision in five space- 
time dimensions, together with the relevant physical quantities. We show that the total radiated 
energy, when two black holes collide from rest at infinity, is approximately (0.089 ± 0.006)% of the 
centre of mass energy, slightly larger than the 0.055% obtained in the four dimensional case, and 
that the ringdown signal at late time is in very good agreement with perturbative calculations. 

PACS numbers: 04.25. dg, 04.50. Gh 



I. INTRODUCTION 

Black objects in higher dimensional space-times have a 
remarkably richer structure than their four dimensional 
counterparts. They appear in a variety of configura- 
tions (e.g. black holes, black branes, black rings, black 
Saturns), and display complex stability phase diagrams 
(see for an overview). They might also play a key 
role in high energy physics: High energy physics scenar- 
ios, such as the gauge-gravity duality 0] or TeV gravity 
models , suggest that dynamical processes involving 
higher dimensional black holes (BHs) may be relevant for 
understanding the physics under experimental scrutiny 
at particle colliders, such as the Large Hadron Collider 
(LHC) or the Relativistic Heavy Ion CoUider (RHIC). 
These processes can be quite violent and highly nonlin- 
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ear, as in the case of BH collisions. Numerical relativ- 
ity, which solves Einstein's equations on supercomput- 
ers, is therefore the only available tool for high-precision 
studies of such BH systems. Fortunately, the field of 
numerical relativity has matured considerably over the 
last five years (see Refs. d, for reviews), and its 
techniques can now be extended to a much wider class 
of space-times. Space-times of generic dimensionality 
or with more general asymptotics, feature most promi- 
nently among such generalizations. Recent applications 
to the study of higher dimensional BH instabilities may 
be found in Refs. [ill, IHIj using the formalism devel- 
oped in Ref. [l3| ; an application to the study of AdS-\\k.e 
asymptotics may be found in Ref. [l^ |. Further appli- 
cations of numerical relativity to more general types of 
space-times have been discussed in Ref. [1| , hereafter de- 
noted as Paper I, wherein we have started a long-term 
effort to evolve BH space-times in higher dimensions nu- 
merically, and developed a framework to perform numer- 
ical simulations of D dimensional space-times with an 
SO{D-2) isometry group (for L> > 5) or S0{D-2,) (for 
D > 6). 

One scenario in which BH collisions play a very well 
defined role is that of TeV-scale gravity, i.e. scenarios 
in which the fundamental Planck scale is of the order 
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of the TeV. The beginning of the scientific runs at the 
LHC makes accurate theoretical modeUing of the experi- 
mental signatures of this scenario for LHC collisions very 
timely. In this scenario, for centre of mass energies well 
above the TeV threshold, recall that LHC collisions will 
reach 14 TeV, parton-parton collision will be dominated 
by the gravitational interaction, and should be well de- 
scribed by any classical gravitational objects with the 
same gravitational energy. For modelling simplicit y, i t 
is convenient to choose these objects to be BHs p^l - l27l |. 
Due to the dominance of the gravitational interaction, 
we may further neglect the electric charge of the holes; 
charge dependant effects should give sub-leading correc- 
tions to the relevant observables. For sufficiently small 
impact parameters, these trans-Planckian collisions are 
expected to form a BH as follows from Thome's 

hoop conjecture [28j which recently has received sup por t 
from the numerical work of Choptuik and Pretorius [29[ . 
Therefore, there is substantial evidence that modelling 
the individual partons as BHs is not biasing the final re- 
sult of the parton scattering process towards BH forma- 
tion. Indeed, this is the idea that, above the fundamental 
Planck scale, matter does not matter; it only matters the 
gravitational energy that each parton carries. 

After formation, the BH should then decay via Hawk- 
ing evaporation. In order to filter experimental data, 
this process has been modelled by dedicated Monte Carlo 
event generators, such as Truenoir, Catfish, Charyb- 
Dis2 or Blackmax [itI. l30l - [33j . The latter two are being 
used by the ATLAS experiment at the LHC. These gen- 
erators clearly exhibit the very distinct experimental sig- 
natures of BH evaporation, including higher multiplicity 
of jets and larger transverse momentum than those pro- 
duced by any standard model process [s^l- The event 
generators also model the BH production phase from the 
parton-parton scattering, for which they need as an in- 
put the threshold impact parameter for BH formation 
and the energy lost in gravitational radiation during the 
parton-parton collision. At the moment, the best es- 
timates for these quantities are based on trapped sur- 
face methods [s^. Accurate results can, however, be ob- 
tained from full-blown numerical simulations, as has al- 
ready been seen in four dimensions [36l - [38| . Such results 
will be instrumental in a more accurate phenomenolog- 
ical modelling of BH production/evaporation in particle 
colliders. Observe that, even if no evidence for BH for- 
mation/evaporation is found at the LHC, such accurate 
modelling will matter for setting precise lower bounds on 
the fundamental Planck scale. 

In this paper we present the first fully nonlinear treat- 
ment of a head-on collision of BHs in a higher dimen- 
sional space-time, together with an analysis of the rele- 
vant physical quantities. This is achieved by solving the 
corresponding Einstein equations numerically, for which 
we use the formalism and code reported in Paper I. Once 
given the numerically constructed space-time, one is still 
left with the question of extracting physically meaning- 
ful quantities, such as the energy and linear and angular 



momentum carried away by the gravitational radiation. 

In four space-time dimensions, two distinct formalisms 
have been dev elop ed to extract the physical information 
(see e.g. Ref. [33] for a review of both formalisms). One 
is based on the Regge- Wheeler- Zerilli perturbation the- 
ory for the Schwarzschild BH |4]]|: the other is based 
on the Newman- Penrose formalism [42| . which was used 
by Teukolsky to study perturbations of algebraically spe- 
cial space-times [1^, a class which includes the Kerr 
BH. A higher dimensional generalization of the Newman- 
Penrosc formalism has been developed in Refs. [44l - l47| . 
Unfortunately, the condition of being algebraically spe- 
cial does not seem to be as powerful for the study of exact 
solutions or their perturbations in higher dimensions as 
it was in four dimensions. For instance, the Goldberg- 
Sachs theorem is not valid any longer in higher dimen- 
sions 0, For the study we present herein, however, 
it suffices to use the higher dimensional generalisation 
of the Regge- Whccler-Zerilli formalism, because the final 
result of a head-on collision of two D dimensional, non- 
spinning BHs approaches, at late times, a D dimensional 
Schwarzschild, i.e. Tangherlini [4^ BH. Fortunately, the 
perturbation theory of the latter BH has been fully devel- 
oped, in arbitrary dimensions, by Kodama and Ishibashi 
[49| . Our remaining task is to obtain the relevant gauge- 
invariant quantities from our numerical data, a procedure 
that we shall describe in detail in this paper. 

After implementing this wave-extraction formalism we 
apply it to study the head-on collision of two BHs in four 
and five dimensional space-times. In four dimensions we 
recover previous results in the literature, and we perform 
a number of tests on the numerical coordinate system, to 
ensure it is appropriate for the wave extraction formal- 
ism. We estimate that around 0.055% of the centre of 
mass energy is radiated when two BHs, at rest at infinity, 
collide. This result is in good agreement with those re- 
ported in the literature [50| . The five dimensional results 
arc entirely new. We show that the kinematics of the 
BHs before the merging follow, to a good precision, the 
Newtonian prediction. We estimate that around 0.089% 
of the centre of mass energy is radiated as gravitational 
waves, when two BHs collide from rest at infinity, and 
present the associated waveforms. We stress that these 
results refer to a fully nonlinear evolution of Einstein's 
field equations. 

This paper is organized as follows. In Section|lI]we ap- 
ply the Kodama-Ishibashi (KI) formalism to the space- 
times considered in Paper I. By assuming that the numer- 
ical space-time is a small deviation from the Tangher- 
lini solution, one is able to relate (see relations (jll C[) ) 
the numerical metric to the KI metric perturbations and 
to compute gauge- invariant quantities using Eqs. (j2.28p . 
(|2.29p . These are then used to construct a master func- 
tion $ (see Eq. (|2.54[) ') from which all relevant infor- 
mation about the radiation can be computed. In Sec- 
tions mil IIVI we present results obtained from the evolu- 
tion of Brill-Lindquist initial data in I? = 4, 5 respec- 
tively, that represents the collision of two equal-mass. 
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nonspinning BHs which are initially at rest. In order 
to calibrate the accuracy of the wave extraction formal- 
ism, we perform a number of tests, including tests on 
the numerical coordinates themselves. We compute the 
time derivative of the master function $, energy fluxes 
and total energy radiated. We give some final remarks in 
Sec. IVland discuss future steps in this research program. 
Two appendices cover technical details on the coordinate 
transformation between the numerical and wave extrac- 
tion frames (Appendix and on the D dimensional har- 
monic expansion of axisymmetric tensors (Appendix [BJ . 



II. GRAVITATIONAL WAVE EXTRACTION IN 
D DIMENSIONAL AXIALLY SYMMETRIC 
SPACE-TIMES 

A. Coordinate frames 

In the approach developed in Paper I, we perform 
a dimensional reduction by isometry on the {D — 4)- 
sphere S^~^, in such a way that the D dimensional 
vacuum Einstein equations are rewritten as an effective 
3-1-1 dimensional time evolution problem with source 
terms that involve a scalar field. The evolution equa- 
tions are expressed in the Baumgarte-Shapiro-Shibata- 
Nakamura (BSSN) formulation [5ll [52|. and numerically 
implemented using a modification of the Lean code [50[ . 

In Paper I we considered different generalizations of 
"axial symmetries" to higher dimensions: either D > 5 
dimensional space-times with SO{D — 2) isometry group, 
OT D > 6 dimensional space-times with SO{D — 3) isom- 
etry group. In this work we only study the former case, 
which allows us to model head-on collisions of nonspin- 
ning BHs; we dub hereafter these space-times as axially 
symmetric. Although the corresponding symmetry man- 
ifold is the (D — 3)-sphere S^~^, the quotient manifold in 
our dimensional reduction is its submanifold S^~'^. The 
coordinate frame in which the numerical simulations are 
performed is 



(x^ 0\ . . . , 0^-4) = (<, x,y,z,cf,\..., ,^^-4) , (2.1) 

where the angles 0^ , . . . , 0^"^ describe the quotient man- 
ifold S^~'^ and do not appear explicitly in the simula- 
tions. Here, z is the symmetry axis, i.e. the collision line. 

In the frame (j2.ip , the space-time metric has the form 
(cf. Eqs. (2.14) and (2.21) of Paper I) 

ds^ =gf,^{x°')dx''dx'' + X{x'')dnD-4 



- a^dt^ + 7y(da;* + P'dt){dx^ + (3^ dt) 

+ \{x'')dnD~4 , 



(2.2) 



where 



(i,a;*), A(x^) is a scalar field and a,/3* are 



the lapse function and the shift vector, respectively. It is 
worth noting that, although in D = 4 a general axially 
symmetric space-time has nonvanishing mixed compo- 
nents of the metric (like 5*0), in > 5 such components 



vanish in an appropriate coordinate frame, as we have 
shown in Paper I. 

With an appropriate transformation of the four dimen- 
sional coordinates x'^, the residual symmetry left after 
the dimensional reduction on S^~^ can be made mani- 
fest: x^' {xf',9) (^ = 0,1,2), 

gf,^{x°')dx''dx'' = gp,i,{x")dx''dx'' + gee{x")de^ (2.3) 



and 



A(x^) =sin2 9gee{x''), 



(2.4) 



so that Eq. (|2.2p takes the form ds^ = g^i^dx^'dx'^ + 
ggsdriu^s, as discussed in Paper I. 

To extract the gravitational waves far away from the 
symmetry axis we employ the KI formalism [49[, which 
generalizes the Regge-Wheeler-Zerilli [40,|4l| approach to 
higher dimensions. We require that the space-time, far 
away from the BHs, is approximately spherically sym- 
metric. Note, that spherical symmetry in D dimensions 
means symmetry with respect to rotations on 5*^^^; this 
is an approximate symmetry which holds asymptotically 
far away from the axis and which is manifest in the co- 
ordinate frame: 



(x^0,0,0^...,0^-^) = (^,r, 



(2.5) 



Note that = t,r and that we have introduced polar- 
like coordinates 9,9 S [0,7r] to "build up" the manifold 
S^^^ in the background, together with a radial spher- 
ical coordinate r, which is the areal coordinate in the 
background. 

The coordinate frame (|2.5p is defined in such a way 
that the metric can be expressed as a stationary back- 
ground (ds*^°^)^ (i.e., the Tangherlini metric) plus a per- 
turbation (ds^^-')^ which decays faster than l/r^~^ for 
large r: 



^gl?df + gi°^dr^ + 
=g[t''dt^ + gi^Jdr^ + {dP + sin^ 9dn 



D-3 



D-3 



D-3 



dr^ 



+ [d9^ + sin^ 6 {d9^ + sin^ 9dnD-4)] , 

(2.6) 

{ds'^^yf ^habdx''dx'' + h^gdx''d9 + hggdS^ 

+ hggdnD-3- (2.7) 

Here, the Schwarzschild radius rs replaces the parameter 
/i used in Paper I and is related to the Arnowitt-Deser- 
Misner mass M by 



„D-3 



IGttM 



{D - 2)Ad-2 ' 



(2.8) 
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where Ad-2 is the area of the {D — 2)-sphere (see 
Eq. ((B20)) ). For instance, = 2M in D = 4 and 
rs = ^/SM/JS^ in L> = 5. 

When we define the coordinate frame (|2.5p . we also 
require that the coordinate 9 in this frame coincides with 
the coordinate 9 appearing in Eq. (j2.3p . With this choice, 
the axial symmetry of the space-time implies that 



0. 



(2.9) 



as in Eq. (PT7)) . and A = sin^ 9ggg, i.e. Eq. (j^ . 

The transformation from the coordinates a;^ = 
(t, X, y, z) in which the numerical simulation is imple- 
mented to the coordinates (x", 9, 9) = (f, r, ^, 6*) in which 
the wave extraction is performed is given by 



X = R sin 9 cos 9 , 
y = R sin sin , 
z = R cos ^ , 



(2.10) 
(2.11) 
(2.12) 



where R = x"^ + y'^ + z-^ and by the reparametrization 
of the radial coordinate 



R = R{r) 



(2.13) 



We note that Eqs. (|2.12p . (|2.13p correctly transform 
the three-metric 7,,, describing the initial data 

-l.jdx'dx' = ipT^ {dR^ + R^{d9^ + sin^ 9d9^)) , (2.14) 

where, as discussed in Paper I, we choose Brill-Lindquist 
initial data. 



?A = 1 + 



' S,2 



4r 



n-3 



4r. 



n-3 



(2.15) 



in order to simulate a head-on collision of BHs starting 
from rest. Here, rs^i and denote the Schwarzschild 
radius and the initial position of the i-ih BH, respec- 
tively. Far away from the axis the conformal factor is 
given by '0 — >■ 1 -|- perturbations. Therefore, the split- 
ting of the metric into a Tanghcrlini background plus 
a perturbation, Eqs. (j2.6|) - (|2.7p . can be recovered on the 
initial time-slice, if we define the reparametrization (|2.13p 
appropriately. 

Our guess is that the transformation (|2.12p . (|2.13p 
yields the "Tangherlini-|-perturbation" splitting (|2.6p . 
p.7p during the entire evolution of the system. This 
statement can be checked numerically by verifying the 
following relations (see Appendix [B|) : 



d9 sin 



D-3 , 



d9gu{9,9)-gf. 



(0) 



Qtr = 



d9sm^-^9 I d9g, 



1 



d9sm^~^9 



RR( 



= 0. 



0, 
(2.16) 



(2.17) 



q(") = 
(2.18) 



where K^'^ = d9{sm9)^ ^, together with the axisym- 
metry conditions (|2.4|) . (|2.9p . As we will discuss in Sec- 
tion Hill Eqs. (1231), (1221), (|2l6| - ([24^ are indeed sat- 
isfied with high accuracy throughout the numerical evo- 
lution. The preservation of the above identities during 
the numerical evolution justifies also the identification of 
the time coordinate in the numerical and wave extraction 
frames, and our use of the KI formalism. 

Finally, Eqs. (123), (1221), (123) yield the 3 + 1 splitting 

ds^ =(ds(°))2 + (ds(i))^ 

=gfipdxf^dx^ + (r^ sin^ 9 + hge)dnD-3 
=gp,i^dx'^dx^ + (r^ sin^ + hgg) 

X {d9^ + sin^ edno-i) 
= - a^dt^ + 7y [dx' + P^dt) 

X {dx^ + P^dt) + \dVLD-i , (2.19) 

where x'^ = {t, r, 9). With the 3-1-1 splitting, the axisym- 
metry conditions (|2.4p . (|2.9p take the form 



A = jee sin 9 , jug 



(2.20) 



The variable r can be determined from the angular com- 
ponents of the metric (|2.19p . by averaging out hgg, hgg 
(see Appendix |B| ; its explicit expression is given by 



{r{R)f = 



+{D-3)jgg{sm9f~^ 



d9 [7ee(sin6l) 



D-3 



(2.21) 



As we will discuss in Section UlTl we find that the areal 
radius r is very close to R. 



B. Harmonic expansion 

In the KI formalism [i^ (see also [11]), the background 
space-time has the form (|2.6p 



(ds("))2 ^g'^^'gdx'^dx" ^ g'-^^dx'^dx" + r'dnD-2 



JO) 



gffjdx^dx'' + r'-fi-jd(t)'d(ty> , 



(2.22) 



i.e. the Tanghcrlini metric, where the x^ coordinates 
refer to the full space-time. The space-time perturba- 
tions can be decomposed into spherical harmonics on the 
[D - 2)-sphere 5^"^. They are functions of the D - 2 
angles cjf' = {9,9,(j)-^, . . . , We denote the metric of 

S^~^ by 7ij, and with a subscript .j the covariant deriva- 
tive with respect to this metric. Finally, we denote the 
covariant derivative with respect to the metric g^ij with 
a subscript \a- 

As discussed in [4^ , there are three types of spherical 
harmonics: 

• The scalar harmonics 5(0*), which are solutions of 
as = i^S^ = -k^S , (2.23) 
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with k"^ = l{l + D -3), I = 0,1,2,... . The scalar 
harmonics S depend on the integer I and on other 
indices; we leave such dependence implicit. We also 
define 



1 



fc2 :y £)_2 



Ti-jS. 



(2.24) 



Observe that ^^■'S-fj = 0. 

Each harmonic mode of the metric perturbation 
SgMN ~ can be decomposed as 



Sgab — hab — fabS , 



(2.25) 
(2.26) 

% = = 2r2 {HlJi^S + HtS^^ , (2.27) 

where /ah, fa, Hl, Ht are functions of x"" = {t, r). 
Note, that in each of these expressions there is a 
sum over the indices of the harmonic. 

For ^ > 1, the metric perturbations can be ex- 
pressed in terms of the following gauge-invariant 
variables [531 



F = Hl + j^Ht + i^ar-l" , 
D — 2 r 

Fab = fab + Xa\b + , 

where we have defined 

r / r 
Xa = - { fa + yHria 



(2.28) 
(2.29) 



k V^" k 

• The vector harmonics V:;(0*), solutions of 

Y'V^,^ = -4Vfc , (2.30) 

with k'^ = l{l + D - 3) - I, I = 1,2, ... . These 
harmonics satisfy the relation 



V.*, = . 



(2.31) 



The harmonic expansion of the corresponding met- 
ric perturbations is given by Eqs. (|2.26p - (|2.27p . 
with Si replaced by V^, replaced by 



(2.32) 



and Hl = 0. 
• The tensor harmonics Tjj (</>*), solutions of 

7*"'%g:ij = —k^Tfs , (2.33) 



with fcf, ?(/ -f D - 3) - 2, ; = 1, 2, , 
harmonics satisfy, 



0. 



These 



(2.34) 



In the Z? = 4 case they vanish. The harmonic ex- 
pansion of the corresponding metric perturbations 
is given by (|2.27p . with S^-j replaced by Tij and 

i?L = 0. 



C. Implementation of axisymmetry 

In an axially symmetric space-time, the metric pertur- 
bations are symmetric with respect to S^^^ . Therefore, 
the harmonics in the expansion of Iimn depend only on 
the angle (which does not belong to S^~^). Further- 
more, since there are no off-diagonal terms in the metric 
(cf. Paper I), the only nonvanishing g^^ components are 
5^^; the only components gij are either proportional to 
Jlj, or all vanishing but ggg. This implies that only scalar 
spherical harmonics can appear in the expansion of the 
metric perturbations. Indeed, if 

V* = (V^0,...,0), V- = V\9), (2.35) 
then Eq. ([23T|) gives 

"1^:1 = "^1 = ^ = ^ = . (2.36) 



Similarly, from Eq. (|2.34p we obtain 7ij = 0. 

The scalar harmonics, solutions of Eq. (|2.23p and which 
depend only on the coordinate 9, arc given by the Gcgen- 
bauer polynomials C^^^^^^^, as discussed in Refs. [s^ - 
[56j : writing explicitly the index /, they take the form 

Si{e) = {K^'^)-^l^c\°'^^'^{cose) , (2.37) 
where the normalization K^^ is chosen such that 



dn^-^Si^gSi,^g = Su'k^ , 
(2.38) 

and k^ ^ l{l + D — 3) (see Appendix |B]). By computing 
Si J, Saj from Eqs. (p:^ (using Eq. (P?^ ') we find 



D — 3 
k^{D^2)^'' 
sir? 9 



(2.39) 
(2.40) 



fc2(L»-2) 
where we have defined 

m{9) = Si cot eSi-g. (2.41) 
Therefore, the metric perturbations are given by 

hab=fibSi{9), (2.42) 
Ks =rf'aSi{0)e = '\rflSi{9)-s , (2.43) 
hee =2r\HiSi{9) + HlrSi{9)ss) 

=2r^ (niSiie) + H}rj^^—^Wii9)] , (2.44) 



k^iD-2) 
hee =2r2(if[ sin^ 0Si{9) + Hi^Si{9)ge) 



2r^ sm^9[ HlSi{9) - iJ^ 



B{D-2) 



m{9) 



(2.45) 
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The quantities fab, fa, Hl, Ht are (see Appendix IB|) 



^ J\eisiner-^K,Cr'^/\ (2.46) 



fa{t,r) 



1 



Ad-3 



y/lil + D-3)rVK^ 



(2.47) 



HLit,r) = 



1 Ad-3 
2{D - 2)r2 ^ 



d9 {sin 9) 



D-3 



ha 



D-3 



ho 



sm 



HT{t,r) 



Ad-3 



2r2(fc2 -15 + 2) VX'^ JO 



c(^-^)/^(cos0"), 

(2.48) 



hao 



Wlie). 



(2.49) 



where hab = hab{t,r,9), h^g = hae{t,r,0), 
^eS = hgg{t,r,9), hgg = hgg{t,r,9) and c\^~^^^'^ = 
c/^-^'/^(cos0-). 

In terms of these quantities, using Eqs. (|2.28p . (|2.29p , 
we get the gauge- invariant quantities F, Fab- 

As we have discussed above, this approach has been de- 
veloped for D > 4, since in D = 4 the off-diagonal terms 
gtip, drip are not vanishing in general axially symmetric 
space-times. However, we can extend our framework to 
Z? = 4 if we restrict ourselves to axially symmetric space- 
times with gtrj, = grrf, = 0. In this way, we can test our 
formalism by comparing our results to the existing liter- 
ature. For instance, we note that in Z? = 4 the perturba- 
tion functions are related to the expressions in Ref. [st} . 
with the identifications 



/afc — -^0, Hi, H2, 
^j:-fa ~ ho, hi , 



2Hr 



2Ht 



= K . 



(2.50) 
(2.51) 

(2.52) 
(2.53) 



We also remark that in the transverse-traceless gauge, 
only Ht is nonvanishing, but in a generic gauge (like the 
one used in the numerical simulations) all these quantities 
are in principle nonvanishing. 



D. Extracting gravitational waves at infinity 

In the KI framework, the emitted gravitational waves 
are described by the master function To compute 
$ in terms of the gauge-invariant quantities F , Fab one 
should perform a Fourier transform or a time integration 
(see [49|). This can be avoided if we compute directly 



4>,t, given b£| 

= (i?-2)r(^-4)/2. 



-F\ + 2rFt 



A:2 - D + 2 



(D-2)[D-1) ' 



(2.54) 

where k"^ = /(/ + £>- 3). In the TT-gauge, the gravi- 
tational perturbation is described by Ht, which decays 
as r('^~2)/2 ^ffll}l increasing r, whereas the other pertur- 
bation functions have a faster decay (see [13] )■ In this 
gauge, the asymptotic behaviour of the master function 
is 



(2.55) 



and tends to an oscillating function with constant am- 
plitude as r — > 00. The asymptotic behaviour of $ has 
been checked numerically (cf. Section llll| . 

Writing the index I explicitly, the energy flux in each 

multipole is [H^ 

The total energy emitted in the process is then 



00 „ + oo 



El , dEi 



1=2 



(2.57) 



III. 



HEAD-ON COLLISION FROM REST IN 

D = 4 



The numerical simulations of head-on collisions of 
equal-mass binaries starting from rest have been per- 
formed with the Lean code originall y i ntroduced in 
Ref. [HO], modified along Sec. Ill of Ref. [HI] and adapted 
to higher dimensional space-times in Paper I. The Lean 
code is based on the Cactus computational toolkit 15911 
and uses the Carpet mesh refinement package l60l. |61| , 
the apparent horizon finder AHFinderDirect [62ll63l| and 
the puncture initial data solver of Ref. . Head-on col- 
lisions in four dimensional space-times have been stud- 
ied extensively in the literature and provide valuable op- 
portunities to calibrate the wave extraction formalism. 
These tests are the subject of the remainder of this Sec- 
tion, while we discuss our new results for five dimensional 
space-times in Sec. II VI below. 

In order to test our implementation of the KI formal- 
ism in D = 4, we have simulated head-on collision of an 
equal-mass, non spinning BH binary initially at rest. The 
parameters used in this simulation are shown in Table HI 
This particular system is well understood and enables us 
to compare our results derived from the KI formalism 



^ Note that there is a factor r missing in Eq. (3.15) of Ref. [49 
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TABLE I. Grid structure and initial parameters of the head-on coUisions starting from rest in D = 4 and D — 5. The grid setup 
is given in terms of the "radii" of the individual refinement levels, in units of rg, as well as the resolution near the punctures 
h (see Sec. II E in [sOl for details) . d is the initial coordinate separation of the two punctures and L denotes the proper initial 
separation. 



Run 


D 


Grid Setup 




d/rs 


L/rs 


HD4c 


4 


{(128, 64, 32, 16, 8) x (1, 0.5, 0.25), h = 


rs/80} 


5.257 


7.154 


HD4„ 


4 


{(128, 64, 32, 16, 8) x (1, 0.5, 0.25), h = 


r-s/88} 


5.257 


7.154 


HD4/ 


4 


{(128, 64, 32, 16, 8) x (1, 0.5, 0.25), h = 


r-s/96} 


5.257 


7.154 


HD5a 


5 


{(256, 128, 64, 32, 16, 8, 4) x (0.5, 0.25), h 


= rs/84} 


1.57 


1.42 


HD5b 


5 


{(256, 128, 64, 32, 16, 8, 4) x (0.5, 0.25), h 


= rs/84} 


1.99 


1.87 


HD5c 


5 


{(256, 128, 64, 32, 16, 8, 4) x (1, 0.5), h = 


= rs/84} 


2.51 


2.41 


HD5d 


5 


{(256, 128, 64, 32, 16, 8, 4) x (1, 0.5), h = 


= rs/84} 


3.17 


3.09 


HD5ec 


5 


{(256, 128, 64, 32, 16, 8) x (2, 1, 0.5), h = 


= rs/60} 


6.37 


6.33 


HD5e™ 


5 


{(256, 128, 64, 32, 16, 8) x (2, 1, 0.5), h = 


= rs/72} 


6.37 


6.33 


HD5e/ 


5 


{(256, 128, 64, 32, 16, 8) x (2, 1, 0.5), h = 


= rs/84} 


6.37 


6.33 


HD5f 


5 


{(256, 128, 64, 32, 16, 8) x (2, 1, 0.5), h = 


= rs/84} 


10.37 


10.35 



against those obtained using both, the Rcggc-Wheeler- 
Zcrilli wave extraction and the Newman-Pcnrose frame- 
work; cf. 57 1 and [s^ for corresponding literature stud- 



ies. 

In order to perform these tests, we need to relate our 
master function $ of Scc. lIIDl to the variables used in tra- 
ditional four dimensional studies. Specifically, a straight- 
forward calculation shows that the Zerilli wavefunction $ 
adopted in Ref. [s^l for Z = 2 multipoles and the outgoing 
Weyl scalar ^'4 used in [s^l can be expressed in terms of 
$ according to 



$ 6$, 
r*4 = ^*,tt 



(3.1) 

(3.2) 



Note that the imaginary part of vanishes in the case 
of a head on collision, due to symmetry. 

The resolution in h — rg/9Q for all results reported in 
this section except for the convergence study in Sec. lIIICl 
which also uses the lower resolutions he = "rs/SO and 
hm = ''s/SsQ Gravitational waves have been extracted 
at three different coordinate radii R (cf. Eq. I2.13P , which 
we denote by i?cx = 30 , 40 rs , 50 rg. 



A. Tests on the numerical coordinates 

The procedure described in Section HIl assumes that the 
numerical space-time consists of a small deviation from 
the Schwarzschild-Tangherlini metric. In order to ensure 



that the gravitational waves are extracted in an appro- 
priate coordinate system we perform a number of checks. 
We first test the relations ((2T6l) . ([2T7)) and ((2T8)) . In 
Fig. [l]we show Qu, i.e., the difference between the numer- 
ical gtt , averaged over the extraction sphere and the cor- 
responding component of the assumed background met- 
ric. Here we evaluate the background metric by assum- 
ing, as a first approximation, that the Schwarzschild ra- 
dius of the BH is rs = rs^i + rs,2- 

The deviation of the full 4-metric from the 
Schwarzschild-Tangherlini background decreases as the 
extraction radius increases. Indeed, a straightfor- 
ward calculation shows that a deviation Srs of the 
Schwarzschild radius from the background value leads to 
Gtt ^ Srg~^ /r^~^, i. e. 6rs/rtor D 4. In the left panel 
of Fig.[l]wc therefore show the deviation Qu rc-scalcd by 
r. We further apply a time shift to account for the differ- 
ent propagation time of the wave to reach the extraction 
radii. As shown in the figure, the deviation from the 
Schwarzschild line element is small and decreases 1/r 
in accordance with our expectation. We also note that a 
deviation Srs represents a monopole perturbation of the 
background which decouples from the quadrupolc wave 
signal at perturbative order, so that its impact on our 
results is further reduced. 

In summary, wc can give an uncertainty estimate for 
the approximation rs ~ rs^i -\-rs^2 f^or the Schwarzschild 
radius of the final BH, which ignores the energy loss 
through gravitational radiation. As demonstrated by the 
left panel of Fig. [U at late times |i?ex/rs Qtt\ ~ 0.01, and, 
since r ~ i?ox (as we discuss below), we obtain the upper 
bound 



^ In order to ensure that our fundamental unit is of physical di- 
mension length for all values of space-time dimension D, we be- 
lieve it convenient to express our results in units of the radius 
„l5-3 D-3 I „D-3n 



rs (given by ' 



'S,i 



+ r 



S,2 



) of the "total" event hori- 



zon as opposed to the total BH mass M commonly used in four 
dimensional numerical relativity. In D = 4, of course, rs = 2M. 



Srs ^ r 



0.01. 



rs 



(3.3) 



This crude analysis sets an upper bound of 1% on the 
fraction of the centre of mass energy radiated as gravita- 
tional waves. We further note that the close agreement 
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FIG. 1. (Color online) Left panel: Qtt calculated from Eq. (|2.16p for D — 4, aX different extraction radii. This quantity has 
been shifted in time to account for the different extraction radii and re-scaled by the corresponding _Rcx. The late time behavior 
is shown in the inset. Right panel: time evolution of areal radius (cf. (|2.2ip ) extracted at the radii i?cx ~ SOrg (black solid 
line), i?ox = 40rs (red dashed line) and Rex ~ 50rs (green dashed-dotted line). 



between gu and its Tangherlini counterpart implies that 
the time coordinate employed in the numerical simula- 
tion and the Tangherlini coordinate time coincide. By 
analysing Qtr and Qrr hi the same manner, we find that 
relations (|2.16p - (|2.18p are satisfied with an accuracy of 
one part in 10^ throughout the evolution, and one part 
in 10'^ at late times, when the space-time consists of a 
single distorted black hole. 

In practice, gravitational waves are extracted on spher- 
ical shells of constant coordinate radius. The significance 
of the areal radius associated with such a coordinate 
sphere in the context of extrapolation of GW signals has 
been studied in detail in Ref. [65[. For our purposes, the 



most important question is to what extent gauge effects 
change the areal radius (|2.2ip of our extraction spheres. 
For this purpose, we show its time evolution in the right 
panel of Fig. [1] for different values of i?ox- The reassur- 
ing result is that the areal radius exceeds its coordinate 
counterpart by about 1 % at i?ox ~ 50 rs and remains 
nearly constant in time. 



B. Waveforms 



cellent agreement between the different extraction meth- 
ods. Next we consider the dependence of the wave signal 
on the extraction radius. In Fig. [3] we show the I = 2 
component of extracted at three different radii and 
shifted in time by i?ox- As is apparent from the figure, 
the wave function shows little variation with i?ox at large 
distances, in agreement with expectations. 

A further test of the wave signal arises from its late- 
time behaviour which is dominated by the BH ringdown 
[66j . an exponentially damped sinusoid of the form e~"^*, 
with uj being a characteristic frequency called quasinor- 
mal mode (QNM) frequency. Using well-known meth- 
ods [66| - |68| , we estimate this frequency to he rs oj ~ 
0.746 ± 0.002 - i (0.176 ± 0.002). This can be com- 
pared with theoretical predictions from a linearized ap- 
proach, yielding rguj = 0.747344- i 0.177925. Finally, 
we consider the numerical convergence of our results. In 
Fig. 01 we plot the differences obtained for ^ t extracted 
at i?ex = 30 rs, using the different resolutions of the three 
models HD4 listed in Table HI The differences thus ob- 
tained are consistent with 4*'' order convergence. This 
implies a discretization error in the / = 2 component of 
$.( of about 4% for the grid resolutions used in this work. 



As a benchmark for our wave extraction, wc compare 
our results obtained with independent wave extraction 
tools; (i) the explicitly four dimensional Zerilli formalism 
and (ii) the Newman-Penrose scalars. For this purpose 
we have evolved model HD4/ and extracted the Zerilli 
function according to the procedure described in (STj (see 
also Eqs. (|2.50p - (|2.53p above) and the Newman Penrose 
scalar ^'4 as summarized in |50[. These are compared 
with the KI wave function ^ t and its time derivative $_tt 
in Fig. [21 Except for a small amount of high frequency 
noise in the junk radiation at t w 25'r5, we observe ex- 



C. Radiated energy 

Once the KI function $,f is known, the energy flux 
can be computed from Eq. (j2.56p . For comparison, we 
have also determined the flux from the outgoing New- 
man Penrose scalar 'i'4 according to Eq. (22) in Ref. [13]. 
The flux and energy radiated in the I = 2 multipole, ob- 
tained with the two methods at i?cx — 50 rs is shown in 
Fig. 151 and demonstrates agreement within the numerical 
uncertainties of about 4 % for either result. We obtain 
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FIG. 2. (Color online) Left panel: Time derivatives of the I = 2 modes of the KI function $ (black solid line), and of the 
Zerilli function <1> (red dashed line) extracted for model HD4/ at 7?ex = 30rs- The KI function has been re-scaled by a constant 
factor (cf. (|3.1|l 'l which accounts for the different normalizations of both formulations. Right panel: comparison of the second 
time derivative <I>,tt with the outgoing Newman-Penrose scalar »I'4 for the same model. The KI wavefunction has been re-scaled 
according to Eq. (|3.2|l . 
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FIG. 3. (Color online) The I = 2 component of the KI wave 
function $,t extracted at the radii i?ex = 30rs (black solid 
line), i?cx = 40rs (red dashed line) and _Rcx = 50rs (green 
dashed-dotted line). They have been shifted in time by the 
corresponding Tiex. 



FIG. 4. (Color online) Convergence analysis of the I = 2 
component of extracted at _Rex = 30 rg. We plot the dif- 
ferences between the low and medium resolution (black solid 
line) and medium and high resolution (red dashed line) run. 



The latter is re-scaled by the factor Q 
order convergence [Til ]. 



1.58 expected for 4 



an integrated energy of 5.5 x 10~^ M and 5.3 x 10~^ M, 
respectively, for the gravitational wave energy radiated 
in ? = 2, where M denotes the centre of mass energy. 

The energy in the I = 2 mode is known to contain 
more than 99% of the total radiated energy [5^. Our 
analysis is compatible with this finding; while the energy 
in the I = 3 mode is zero by symmetry, our result for the 
energy in the / = 4 mode obtained from the KI master 
function is three orders of magnitude smaller than that 
of the I = 2 contribution. 



IV. HEAD-ON COLLISION FROM REST IN 

D = 5 

Having tested the wave extraction formalism in four di- 
mensions, we now turn our attention to the new results 
obtained for head-on collisions of BHs in five dimensional 
space-times. As before, we consider nonspinning BH bi- 
naries initially at rest with coordinate separation d. Note 
that in five space-time dimensions the Schwarzschild ra- 
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FIG. 5. (Color online) Energy flux (left panel) and radiated energy (right panel) for the I = 2 mode extracted at Tic: 
from the Kl wave function (black solid curve) and the Newman Penrose scalar >I'4 (red dashed curve). 



50rs 



dius is related to the ADM mass M via Eq. ()2.8 



8M 



(4.1) 



We therefore define the "total" Schwarzschild radius rs 



such that 



's,i 



' S,2- 



By using this definition, rs 



has physical dimension of length and provides a suitable 
unit for measuring both, results and grid setup. 

As summarized in Table U wc consider a sequence of 
BH binaries with initial coordinate separation ranging 
from d = 3.17 rs to d = 10.37 r^. The table further lists 
the proper separation L along the line of sight between 
the holes and the grid configurations used for the indi- 
vidual simulations. 



Qtt are well below 1 part in 10^ at i?ox = 60 rs- Fur- 
thermore, we can estimate the crudeness of the approx- 
imation r| — I + rg 2 for the Schwarzschild radius of 
the final BH: as shown in the right panel of Fig. |6l at 
late times \Rlx/''"'sStt\ ^ 0.01; this value gives an upper 
bound on the radiated energy. 

For the third test, we recall that our higher dimensional 
implementation does not employ the full isometry group 
of the sphere in D = 5 dimensions and axial symmetry 
manifests itself instead in the conditions (|2.20[) on the 
metric components and the scalar field. We find these 
conditions to be satisfied within 1 part in 10^ and 1 part 
in 10^^, respectively, in our numerical simulations which 
thus represent axially symmetric configurations with high 
precision. 



Tests on the numerical coordinates 



B. Newtonian collision time 



In order to verify the assumptions underlying our for- 
malism, we have analysed the coordinate system in anal- 
ogy to Sec. nil Al First, we have evaluated the averaged 
areal radius on extraction spheres of constant coordinate 
radius. 

The result shown in the left panel of Fig. IH] demon- 
strates that the coordinate and areal radius agree within 
about 1 part in 10"* for i?ox > 40 rs- The Tangherlini 
coordinate r equals by construction the areal radius and 
our approximation of setting r « R^x in the wave extrac- 
tion zone is satisfied with high precision. 

Second, we evaluate the deviation of the metric com- 
ponents according to Eqs. (|2.16p - (|2.18p . From the dis- 
cussion in Sec. IIII Al we expect Gtt ~ '^^/'^l in D = 5. 
Our results in the right panel of Fig. [6] confirm this ex- 
pectation and demonstrate that our space-time is indeed 
perturbatively close to that of a Tangherlini metric at 
sufficient distances from the black holes; deviations in 



An estimate of the time at which the holes "collide" , 
can be obtained by considering a Newtonian approxima- 
tion to the kinematics of two point particles in Z3 = 5. 
In the weak-field regime, Einstein's equations reduce to 
"Newton's law" a = -VS(x), with ft-oo = -2S(x) = 
rg~^/2r^~^. The Newtonian time it takes for two point- 
masses (with Schwarzschild parameters rs^i and ts^t) to 
collide from rest with initial distance L in D dimensions 
is then given by 



^free-fall 

rs 



(4.2) 



where ^ = + n ^ and 



S,2 




D-3 ' 



r(i + i5^) 



(4.3) 
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FIG. 6. (Color online) Left panel: Time evolution of the areal radius r in units of the extraction radius averaged over coordinate 
spheres at -Rex = 20 rs (black solid), 40 rs (red dashed) and 60 rs (green dash-dotted curve). Right panel: Deviation of the 
metric component Ri^/r%Qu calculated from Eq. (|2.16|l at the same extraction radii and shifted in time to account for differences 
in the propagation time of the wave signal. 



For D — A, one recovers the standard result tfrGo-faii = 
TTy^L^/rgrs , whereas for D = 5 we get 

tfrcc-faii = {L/rsf rs . (4.4) 

In general relativity, BH trajectories and merger times 
are intrinsically observer dependent quantities. For our 
comparison with Newtonian estimates we have chosen 
relativistic trajectories as viewed by observers adapted 
to the numerical coordinate system. While the lack of 
fundamentally gauge invariant analogues in general rela- 
tivity prevents us from deriving rigorous conclusions, we 
believe such a comparison to serve the intuitive interpre- 
tation of results obtained within the "moving puncture" 
gauge. Bearing in mind these caveats, we plot in Fig. [7] 
the analytical estimate of the Newtonian time of colli- 
sion, together with the numerically computed time of 
formation of a common apparent horizon. Also shown 
in Fig. [7] is the time at which the separation between the 
individual hole's puncture trajectory decreases below the 
Schwarzschild parameter rs- The remarkable agreement 
provides yet another example of how well numerically 
successful gauge conditions appear to be adapted to the 
black hole kinematics. It is beyond the scope of this pa- 
per to investigate whether this is coincidental or whether 
such agreement is necessary or at least helpful for gauge 
conditions to ensure numerical stability. SufBce it to say 
at this stage that similar conclusions were reached by An- 
ninos et al. (69j and Lovelace et al. [70| in similar four 
dimensional scenarios. 



C. Waveforms 

We now discuss in detail the gravitational wave signal 
generated by the head-on collision of two BHs in five di- 
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FIG. 7. (Color online) Estimates for the time it takes for 
two equal-mass BHs to collide in -D = 5. The first estimate 
is given by the time tcAH elapsed until a single common ap- 
parent horizon engulfs both BHs (diamonds), the second es- 
timate is obtained by using the trajectory of the BHs, i.e., 
the time ttraj at which their separation has decreased below 
the Schwarzschild radius (circles). Finally, these numerical 
results are compared against a simple Newtonian estimate, 
given by Eq. (|4.4|l (blue solid line). 

mensions. For this purpose, we plot in Fig. [8] the / ~ 2 
multipole of the KI function $ i for model HDSe/ ob- 
tained at different extraction radii. Qualitatively, the sig- 
nal looks similar to that shown in the left panel of Fig. [5] 
for D = 4. A small spurious wavepulse due to the ini- 
tial data construction is visible at At w 0. This so-called 
"junk radiation" increases in magnitude if the simulation 
starts with smaller initial separation of the holes. We 
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FIG. 8. (Color online) Left panel: The I = 2 component of the KI waveform for model HD5e/ extracted at radii Ra^/rs = 20, 40 
and 60 and shifted in time by Rcx/rs- Right panel; The I = 2 and I = 4 mode of the KI function for the same simulation, 
extracted at Rcx/rs ~ 60. For clarity, the I = 4 component has been re-scaled by a factor of 100. 



0.003 




FIG. 9. (Color online) Convergence analysis of the I — 2 com- 
ponent of the KI function generated by model HD5e extracted 
at /?ex ~ 60rs- The difference between the medium and high 
resolution waveforms has been amplified by the factors 1.97 
(red dashed line) and 2.33 (green dashed-dotted line) indicat- 
ing third and fourth order convergence. 



return to this issue further below, when we study the de- 
pendence of the gravitational radiation on the initial BH 
separation. The physical part of the waveform is domi- 
nated by the merger signal around At ~ 50 rs, followed 
by the (exponentially damped) ringdown, whereas the in- 
fall of the holes before At = 40 rs does not produce a 
significant amount of gravitational waves. Comparison of 
the waveforms extracted at different radii demonstrates 
excellent agreement, in particular for those extracted at 
Rex = 40 rs and 60 rs- Extrapolation of the radiated 
energy to infinite extraction radius yield a relative error 



of 5 % at i?ex = 60 rs, indicating that such radii are 
adequate for the analysis presented in this work. 

Due to symmetry, no gravitational waves are emitted 
in the I = 3 multipole, so that / = 4 represents the sec- 
ond strongest contribution to the wave signal. As demon- 
strated in the right panel of Fig.[5J however, its amplitude 
is two orders of magnitude below that of the quadrupole. 

A convergence analysis also using the lower resolution 
simulations of models HDSCc and HDSCm is shown in 
Fig. [S] and demonstrates overall convergence of third to 
fourth order, consistent with the numerical implemen- 
tation. From this analysis we obtain a conservative es- 
timate of about 4% for the discretization error in the 
waveform. 

In practice, numerical simulations will always start 
with a finite separation of the two black holes. In order to 
assess how accurately we are thus able to approximate an 
infall from infinity, we have varied the initial separation 
for models HD5a to HD5f as summarized in Table HI For 
small d we observe two effects which make the physical 
interpretation of models HD5a— HD5c difficult. First, the 
amplitude of the spurious initial radiation increases and 
second, the shorter infall time causes an overlap of this 
spurious radiation with the merger signal. As demon- 
strated in Fig. [To] for models HD5e and HD5f, however, 
we can safely neglect the spurious radiation as well as 
the impact of a final initial separation, provided we use 
a sufficiently large initial distance d > 6 of the BH bi- 
nary. Here, we compare the radiation emitted during the 
head-on collision of BHs starting from rest with initial 
separations 6.37 rs and 10.37 rs- The waveforms have 
been shifted in time by the extraction radius i?ex = 60 rs 
and such that the formation of a common apparent hori- 
zon occurs at At = 0. The merger signal starting around 
At = shows excellent agreement for the two configura- 
tions and is not affected by the spurious signal visible for 
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FIG. 10. (Color online) The I = 2 components of the KI 
function as generated by a head-on collision of BHs with initial 
(coordinate) distance d = 6.37 rg (black solid line) and d — 
10.37 rs (red dashed line). The wave functions have been 
shifted in time such that the formation of a common apparent 
horizon corresponds to At = and taking into account the 
time it takes for the waves to propagate up to the extraction 
radius i?ex = 60 rg. 



FIG. 11. (Color online) Energy flux in the I = 2 component 
of the KI wave function $,4, extracted at _Rox = 60 rs, for 
models HDSe/ (black solid line) and HD5f (red dashed line) in 
TableHl The fluxes have been shifted in time by the extraction 
radius i?ex = 60 rs and the time tcAH at which the common 
apparent horizon forms. 



D. Radiated energy 



HDSe at At « -50 rs. 

We conclude this discussion with two aspects of the 
post-merger part of the gravitational radiation, the ring- 
down and the possibility of GW tails. After formation 
of a common horizon, the waveform is dominated by an 
exponentially damped sinusoid, as the merged hole rings 
down into a stationary state. By fitting our results with 
a exponentially damped sinusoid, we obtain a character- 
istic frequency 

rsoj = 0.955 ± 0.005 - i (0.255 ± 0.005) . (4.5) 

This value is in excellent agreement with perturbative 
calculations, which predict a lowest quasinormal fre- 
quency rsuj = 0.9477 - i 0.2561 for / 2 [56, 66, . 

A well known feature in gravitational waveforms gen- 
erated in BH space-times with D = 4 as well as > 4 
are the so-called power-law tails (T^ - lTSj . In odd dimen- 
sional space-times an additional, different kind of late- 
time power tails arises, which does not depend on the 
presence of a BH. These are due to a peculiar behavior of 
the wave-propagation in flat odd dimensional space-times 
because the Green's function has support inside the en- 
tire light-cone [7l|. We have attempted to identify such 
power-law tails in our signal at late-times, by subtracting 
a best-fit ringdown waveform. Unfortunately, we cannot, 
at this stage, report any evidence of such a power-law in 
our results, most likely because the low amplitude tails 
are buried in numerical noise. 



Comparison of Figs. [3] and [TOl for the GW quadrupole 
in Z? = 4 and D = 5 shows a larger wave amplitude in 
the five dimensional case and thus indicates that this case 
may radiate more energy. We now investigate this ques- 
tion quantitatively by calculating the energy flux from 
the KI master function via Eq. ()2.56|) . The fluxes thus 
obtained for the I = 2 multipole of models HD5e/ and 
HD5f in Table U extracted at Rex = 60 rs, are shown 
in Fig. [TT] As in the case of the KI master function in 
Fig. [TUl we see no significant variation of the flux for 
the two different initial separations. The flux reaches a 
maximum value of dE/dt ^ 3.4 x 10^'* rs, and is then 
dominated by the ringdown flux. The energy flux from 
the I = 4 mode is typically four orders of magnitude 
smaller; this is consistent with the factor of 100 differ- 
ence of the corresponding wave multipoles observed in 
Fig. m and the quadratic dependence of the flux on the 
wave amplitude. 

The total integrated energy emitted throughout the 
head-on collision is presented in the left panel of Fig. [T^l 
We find that a fraction of Erad/M = (8.9 ±0.6) x 10"^ of 
the centre of mass energy is emitted in the form of gravi- 
tational radiation. We have verified for these models that 
the amount of energy contained in the spurious radiation 
is about three orders of magnitude smaller than in the 
physical merger signal. 

An independent estimate for the radiated energy can 
be obtained from the apparent horizon area in the 
effective four dimensional space-time by using the spher- 
ical symmetry of the post-merger remnant hole. Energy 
balance then implies that the energy E radiated in the 
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FIG. 12. (Color online) Left panel: Fraction of the centre of mass energy, Erad/M, radiated in the / = 2 mode of the KI function 
shifted in time such that the origin of the time axis corresponds to the formation of a common apparent horizon. Right panel: 
Fraction of the centre of mass energy 1 — Mah /M radiated during the collision, estimated using apparent horizon information. 
The oscillations in this diagnostic quantity have a frequency comparable to the I — 2 quasinormal mode frequency. 



form of GWs is given by 



E 
M 



— = 1 



Mah 
M 



A, 



47rr| 



(4.6) 



where M^h is the apparent horizon mass. The estimate 
E/M is shown in Fig. [Inland reveals a behavior qualita- 
tively similar to a damped sinusoid with constant offset. 
Indeed, by using a least square fit, we obtain a complex 
frequency rg uj ^ 0.97 — i 0.29, again similar to the funda- 
mental I = 2 quasinormal mode frequency (see discussion 
around Eq. (|4.5p ). At late times, 1 — Mah/M asymp- 
totes to 1 - Mah/M - (9.3 ± 0.8) X 10"'' which agrees 
very well with the GW estimate, within the numerical 
uncertainties. 



V. DISCUSSION 

In this paper we have developed a formalism to ex- 
tract gravitational radiation observables from numerical 
simulations of head-on collisions of BHs in D dimen- 
sions. Moreover, we have performed such simulations 
in D = 4, 5. The D = A case serves as a test of our for- 
malism and demonstrates consistency of our results with 
the literature. The D ~ 5 case is entirely new. Besides 
obtaining the corresponding waveforms, we have shown 
that the total energy released in the form of gravitational 
waves is approximately (0.089±0.006)% of the initial cen- 
tre of mass energy of the system, for a head-on collision 
of two BHs starting from rest at very large distances. As 
a comparison, the analogous process in D = 4 releases 
a slightly smaller quantity: (0.055 ± 0.006)%. We sum- 
marize the main results for head-on collisions of two BHs 
starting from rest in four and five space-time dimensions 
in Table HH 



We have further performed a variety of tests of the 
wave extraction formalism. Besides testing the proxim- 
ity of the numerical coordinate system to the Tangher- 
lini background space-time, we have demonstrated good 
agreement between the radiated energy as derived di- 
rectly from the KI master function with the values ob- 
tained from the horizon area of the post-merger remnant 
hole. Finally, the ringdown part of the waveform yields a 
quasinormal mode frequency in excellent agreement with 
predictions from BH perturbation theory. 

The radiative efficiency E'^^'^/M in Table HIl shows that 
head-on collisions starting from rest in five dimensions 
generate about 1.6 times as much GW energy as their 
four dimensional counterparts. It will be very interesting 
to investigate to what extent this observation holds for 
wider classes of BH collisions. We can compare the radi- 
ation efficiency with the upper limit derived by Hawking 
[76| from the requirement that the horizon area must not 
decrease in the collision. This leads to the area bound 
< 1 — 2~-D-2 . Evidently, this bound decreases with 



dimensionality, while in the present computation it in- 
creases when going from _D ~ 4 to £> = 5. As also 
shown in the table, the generation of GWs in head-on 
collisions starting from rest is about 3 orders of mag- 
nitude below this bound. In four dimensions it has al- 
ready been demonstrated that there exist more violent 
processes which release more radiation than the head-on 
collisions considered in this work fSmSg ■ 



In the context of this work, it would be particularly 
interesting to compute the gravitational radiation emit- 
ted when a point particle falls into a higher dimensional 
BH (the four dimensional calculation is done in the clas- 
sic work by Davis et al (Z^l)- This analysis can be done 
by linearizing Einstein's equations. While such an analy- 
sis was done for infall at high energies [13, HH , it has not 
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TABLE II. Main results for head-on collisions in _D = 4 and 5 dimensions. We list the ring down frequency u, the total energy 
radiated in gravitational waves, the upper bound E'^'"^^ on the radiated energy obtained from Hawking's area theorem and the 
fractional energy in the / = 4 multipole relative to the quadrupole radiation. 



D 


rs 


= 2) 




£='"'=^/A/(%) El^l/E'i^i 


4 


0.7473 - 


i 0.1779 


0.055 


29.3 < 10^^ 


5 


0.9477 - 


i0.2561 


0.089 


20.6 < 10"** 



been done for infalls from rest. The four dimensional case 
shows that by scaling the point-particle results properly 
with the reduced mass, one gets surprisingly good agree- 
ment with full nonlinear studies [TSj . An obvious ques- 
tion is whether such an agreement extends to generic 
number of space-time dimensions. Investigations with 
a similar purpose, but using a different technique, were 
carried out in Refs. 

The numbers reported here for the total energy loss in 
gravitational waves should increase significantly in high 
energy collisions, which are the most relevant scenarios 
for the applications described in the Introduction. In- 
deed, in the four dimensional case, it is known that ultra- 
rclativistic head-on collisions of equal mass nonrotating 
BHs release up to 14% of the initial centre of mass energy 
into gravitational radiation [s^ . The analogous number 
in higher dimensions is as yet unknown and will be sub- 
ject of the next stages of our research programme. Pre- 
liminary results by Gal'tsov et al [8ll,[8^ strongly suggest 
enhancement of radiation emission for higher dimensions, 
in agreement with the D = 5 results shown here. Even 
more energy may be released in high energy collisions 
with nonvanishing impact parameter. In 37,|38| it was 
shown that this number can be as large as 35% in Z? = 4. 
The formalism developed in Paper I allows, in principle, 
the study of analogous processes in D > 6. We hope to 
be able to report on these results in the near future. 
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Appendix A: Coordinate transformation 



In order to extract gravitational radiation using the KI 
formalism one has to perform a coordinate transforma- 
tion from Gartesian coordinates, which are used during 
the numerical evolution, to those adapted for wave ex- 
traction. The physical 3- metric jij, the lapse function a 
and the shift vector /3* computed on our Cartesian grid 
are interpolated onto a Gartesian patch. In terms of these 
quantities we compute the 4-metric g^j/ in Cartesian co- 
ordinates according to Eq. 
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- -f.jP^dtdx' + jijdx'-dx^ . (Al) 



Then, we transform the 4-metric in Cartesian coordinates 
into spherical coordinates, defined by Eq. (|2.12p 



a; = i? sin ^ cos 9 , 
y = R sin 6 sin 6 , 
z = R cos 6 , 



(A2) 
(A3) 
(A4) 



where 9,9 G [0, tt] and R = x'^ + y'^ + z'^ . If we denote 
the metric in spherical coordinates by g^^ and define p = 
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^yx^ + y^, the explicit form of the transformation is 



9tR —9tx sin 9 cos 9 + gty sin 6 sin 9 + gtz cos , 


(A5) 


-.ee 


— Si ^00 


9% ^z{gtx cos 6* + gty sin 6*) - pgt^ , 


(A6) 


Si -.00 


— ^ 00' 


sfe = - y9tx + a;gty , 


(A7) 




- -^01 



9rr —9xx sin^ 9 cos^ 6 + 2172:4/ sin^ ^ cos 9 sin 
+ 2g3;2 sin 9 cos 6* cos 9 + sin^ ^ sin^ 9 
+ 25y2 sin 9 sin 6* cos ^ + cos^ 9 , ( A8) 

—^(52:2: sin^cos^ 9 + 2(7j;j, sin ^ cos 9 sm9 

+ sin 9 sin^ 6* + 5^:2 cos 9 cos 61 + gy^ cos ^ sin 6') 
- {xgxz + yffyz + zg:,^) sin^, (A9) 

5fle —{~'y9xx sin ^ cos + xg^y sin ^ cos — ygxy sin ^ sin 9 



(AlO) 



+ x^yy sm 9 sni 6^ — j/gj;^ cos 9 + xgy^ cos ( 

g^g [gxx cos^ 9 + 2pa;.y cos sin 6* + gyy s\t? 9) 

- 2z{xgxz + y9yz) + P^gzz , (All) 

9S0 ^z{-ygxx cos 6* + xg^y cos 6* - yg^y sin 6* 

+ xgyysui9) + p{ygxz- xgyz) , (A12) 

=-R^ sin^ 9{gxx sin^ — 2(7a;y cos sin 9 + cos^ 9) . 

(A13) 

Henceforth, we wiU drop the superscript 5* and use g^^i, 
for the metric in spherical coordinates. 

The areal radius r is related to i? by a reparametriza- 
tion R = R{r), given by Eq. (jB47p . which depends 
on the components ggg, g00 only. As shown in Sec- 
tion mil we find that this reparametrization is nearly 
constant throughout our numerical simulations. There- 
fore, the quantities grr, gtr, 9r0-i 9r9 can be obtained from 
9RR,9tR,gR0,9Re by a simple rescaling: because 



dR 

dr 



1 



(AM) 



we have grr — gRR, and similar relations hold for the 
other components. 



Appendix B: Harmonic expansion of axisymmetric 
tensors in D dimensions 



As discussed in Section HTCl scalar spherical harmonics 
in D dimensions Si{9, 6*, 0^, ... , are solutions of 

Eq. (12:23)) 



aSi=j'^Sia^ = -k'Si 



(Bl) 



with = l{l + D — 3). Axisymmetric scalar spherical 
harmonics are functions of the coordinate 9 only. Si = 
Si{9). Therefore, Eq. ([BT]) becomes 

nSi (9) = Si M + {D-?,) cot 9Si g = -k^Si , (B2) 



(B3) 

VlftSi = sin 9 cos 9Si , (B4) 
r^i^i'^/ ,0 = sin^ 9 sin ^ cos 9Si g , (B5) 

etc. The quantities S^j defined in Eq. (|2.24p are then 
1 



k^{D~2) 
1 



{{D~2)Si,^ + er,jSi) 
diag((i?-3)Wz, 



k^{D-2) 
- sin^ mi , - sin^ 9 sin^ 6* W; , . . 



where 

'Wi{9) = Si gg - cot 9Si g = sin^ 



S, 



sin 9 



(B6) 



(B7) 



Indeed; using Eq. (jB2[) one finds 

k'{D - 2)Si gg ^{D- 2)Si ^gg + k'Si 

^{D -3){Si,gg^ cot 9Si^g), (B8) 
fc2(Z? - 2)Si gg^{D- 2)Si ,00 + fc2 sin^ 9Si 

= SIT? 9{{D - 2) cot 9Si ^g + k^Si) 
= sin^ 9{-Si gg + cot 9Si g) , (B9) 



and therefore 



— 3 

^'''=k^{D-2)^'' 
sin^ 9 
k^{D-2)^' 



(BIO) 
(Bll) 



and likewise for the other components. 

Axisymmetric scalar spherical harmonics, as discussed 
in Sec. Ill CI can be written in terms of Gegenbauer poly- 
nomials (cf. ((2371) '): 



5,(0) = (/f'^)-l/2c(^-3)/2(c 

If we define 



(B12) 



M^; (cos 61) = C;^^-^^^/^(cos6') - cot6'Cj^-"^^/^(cos6') , 

(B13) 

we have 

yVi{9) = (X'^)"1/2m^/^~^'/'(cos^) . (B14) 
We impose the normalization (|2.38p 



dn^-^SiSv = Siv , / dn^-^Si gS^, g = Sii,k 



(B15) 
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Using 

(B16) 



and 



(z + i2_3) (r(i2_3))^r(; + i) 



we have 



where 



D-3 ■ 



A 



^-3 = YTlElf 



(B17) 
(B18) 

(B19) 
(B20) 



is the surface of the {D — 3)-spherc ^. Note that 
Jdn^-^{---) = AD-3jde{sme)^~^---). with the 
definitions (jf^ Su^-^Si j, 



d9{sm9f-^Si{9)Si-{9) ^SwA^U, (B21) 



d9{sm9)^-'Y^SrA'-^ 



(B23) 



d9{sm9)''-'Sio{9)Si,s{9) 

= SwA^l^ . (B22) 
Furthermore, we note that Eqs. (jB2p and (jB7P imply 
Wi + (I? - 2) cot ^5, + k^Si = , 

so that 

D-2 
sin^ 

= .0 + (I? - 2) cot ^W, + (fc2 - + 2)5, , (B24) 
and therefore 



Wi ,s + - 2) cot 9Si M - —^Si -e + k^Si , 



d6'(sin6i)^-^W,Wp 

d^(sin0)^~'''W,sin0 
\ sin 9 



= -{D -2) I d9{sm e)^-^Wi cot 615,, _ 



d9{sin9)^-^Wi,gSi, ^ 



--{k^ -D + 2) / d9{sm9)^^^Si gSi, 
Jo 

--du'Ao^^k'ie -D + 2). 



(B25) 



We thus obtain 

dn^-'^WiWi, = 6wk^{k^ ~D + 2). (B26) 

The perturbations /^,(t,r), fi{t,r), Hi{t,r), Hl^it.r) 
appearing in the expansion of the metric perturbations 

(nfci) 



Kb — fabSl{0) , 

Ke--\rfiSi{9)^, 



(B27) 
(B28) 



h-os = 2t^ ( + gT pp_\) WzW ) , (B29) 



/lee = 2r^ sin^ e W^Si{6) - Wj, 



B{D-2) 



Wi{9) 



(B30) 



arc given by the following integrals, as follows from 
Eqs. (IBT21) . (fell) . (IbTsI) . (ime)) : 



fait,r) 



dn^'-^habSi 



A 



D-3 



ID 



d9{sm9f-^habCl^~^'>/^, (B31) 



1 



dn 



D-2 



kr 

1 Ad-3 

kr y/KiD j„ 



d9{sn^9r-'K,c[y^'\ 



HL{t,r) 



1 



2{D - 2)r 



dn 



D-2 



, D-3, 

haa H TT^ha 



1 



A 



D-3 



2{D - 2)r2 Jo 



hao H . n - h-a 



d9{sm9) 



sin^ei 

D-3 



sm 



C, 



{D-3}/2 
I ■ 



(B32) 
Si 



(B33) 



1 



2r'^{k^ - D + 2) 
1 



dr2 



r>-2 



sm 



Ad-3 



2r^{k^ - D + 2) 



ID 



d9 {sin 9) 



D-3 



sm ( 



Wi. 



(B34) 



where /i, 



ab 



= habit, r,9), h^g = h^g{t,r,9), 



hfjg = hgg{t,r,9), hee = hgg{t,r,9), C\ 
D-3)/: 



(D-3)/2 



cP"^'/^(cos6l) and Wi = Wi{cos9). 



We also note that the background Tangherlini metric 
depends on the I = harmonic only; the integral of its 
components over I > 2 harmonics vanish. Therefore, if 
we decompose the space-time metric (see Appendix 
as g^^i, = gjS + /i^i/ with fi,v = (t, r, 9, 9) and g^^} is 
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the Tanghcrlini background metric, we can compute the 
integrals (jB34p in terms of the metric g^^, 



ftt 

ftr 
frr 



1 A 



D-3 



1 Ad-3 



ID 



d0"(sin0)^-3c(^-^'/' 



d0"(sin0")^-3^(^-^'/^ 



1 A 



D-3 



TT ./kTd 

1 Ad-3 



(B35) 

dOgtrie, 0), 

(B36) 

dOgrrie, 9), 

(B37) 



kr-K ^/K^d 
X / dO{sme)°-'^d-gCf'^^''^ I dOg 



(B38) 



1 A 



D-3 



X / (i^(sin^)^-35eC/^"^'^^ / deg^g{e,e), 



(B39) 



Hr = 



Ht — 



1 



2{D - 2)r^TT y^K^ 

X / de(^ggg{e,e) + 

1 Ad-3 



sin ( 



(B40) 



2(fc2 - D + 2)r^TT y^K^ 
X y" d0(sin^)^-3^, Jd0[gee{ 



sm 



(B41) 



Furthermore, from Eqs. (|Z!28)) and (|B35|) - (|B4ip we 
deduce 

Ft ^dtHL + -^^dtHr + \f{r) (^dtU + ^^t^rffr) , 

(B42) 

F: =f{r) (^frt + ^(aj. + drft) + ^/t 

+ ^{dtHT + rdtdrHT)^ - Idrfir) (/* + IdtHr) 

(B43) 

Conversely, since the perturbations do not depend on the 
I = harmonic, the background metric g^j^ can be ob- 
tained as follows: 



1 



dO sm^--^0 / d0gfti9,0), (B44) 



= 



„(o) 



1 



d0sm"-'^0 I d0gtr{0,0), 



(B45) 



1 



ddsin"-^0 I d0grr{0,0). (B46) 



Finally, to compute the areal radius r we note that 
9eS ^ + hgg and gee = sin^ + hee- Both the pertur- 
bations hgg and hgg contain harmonics of different type 
{Si, Si jj); to extract the background we need the com- 
bination in Eq. (|B40p : 



:- ^—-Krr- [ d0sm"-^0 [ d0 



^ sin^0 



(B47) 



[1] M. Zilhao et ai, Phys. Rev. D81, 084052 (2010), 
[1001.2302]. 

[2] R. Emparan and H. S. Reall, Living Rev. Rel. 11, 6 

(2008), [0801.3471]. 
[3] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998), 

[hep-th/9711200]. 
[4] I. Antoniadis, Phys. Lett. B246, 377 (1990). 
[5] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. 

Lett. B429, 263 (1998), [hep-ph/9803315]. 
[6] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and 

G. R. Dvali, Phys. Lett. B436, 257 (1998), [hep- 

ph/9804398]. 

[7] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 

(1999), [hep-ph/9905221]. 
[8] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 

(1999), [hep-th/9906064]. 
[9] F. Pretorius, 0710.1338. 
[10] I. Hinder, Class. Quant. Grav. 27, 114004 (2010), 

[1001.5161]. 

[11] M. Shibata and H. Yoshino, Phys. Rev. D81, 021501 
(2010), [0912.3606]. 



[12] M. Shibata and H. Yoshino, Phys. Rev. D81, 104035 

(2010), [1004.4970]. 
[13] H. Yoshino and M. Shibata, Phys. Rev. D80, 084025 

(2009), [0907.2760]. 
[14] H. Witek et al, Phys. Rev. D82, 104037 (2010), 

[1004.4633]. 

[15] T. Banks and W. Fischler, hep-th/9906038. 

[16] S. B. Giddings and S. D. Thomas, Phys. Rev. D65, 

056010 (2002), [hep-ph/0106219]. 
[17] S. Dimopoulos and G. L. Landsberg, Phys. Rev. Lett. 

87, 161602 (2001), [hep-ph/0106295]. 
[18] J. L. Feng and A. D. Shapere, Phys. Rev. Lett. 88, 

021303 (2001), [hep-ph/0109106]. 
[19] E.-J. Ahn, M. Ave, M. Cavaglia and A. V. Ohnto, Phys. 

Rev. D68, 043004 (2003), [hep-ph/0306008]. 
[20] E.-J. Ahn, M. CavagUa and A. V. Olinto, Phys. Lett. 

B551, 1 (2003), [hep-th/0201042]. 
[21] A. Chambhn, F. Cooper and G. C. Nayak, Phys. Rev. 

D70, 075018 (2004), [hep-ph/0405054]. 
[22] V. Cardoso, M. C. Espirito Santo, M. Paulos, M. Pimenta 

and B. Tome, Astropart. Phys. 22, 399 (2005), [hep- 



19 



pli/0405056]. 

[23] M. Cavaglia, Int. J. Mod. Phys. A18, 1843 (2003), [hep- 
ph/0210296]. 

[24] P. Kanti, Int. J. Mod. Phys. A19, 4899 (2004), [hep- 
ph/0402168]. 

[25] P. Kanti, Lect. Notes Phys. 769, 387 (2009), [0802.2218]. 
[26] S. N. Solodukhin, Phys. Lett. B533, 153 (2002), [hep- 
ph/0201248]. 

[27] S. D. H. Hsu, Phys. Lett. B555, 92 (2003), [hep- 
ph/0203154]. 

[28] K. S. Thorne, In Magic Without Magic, edited by J. 
Klauder (W.H. Freeman & co., San Francisco, 1972), 231- 
258. 

[29] M. W. Choptuik and F. Pretorius, Phys. Rev. Lett. 104, 
111101 (2010), [0908.1780]. 

[30] M. Cavaglia, R. Godang, L. Cremaldi and D. Sum- 
mers, Comput. Phys. Commun. 177, 506 (2007), [hep- 
ph/0609001]. 

[31] J. A. Frost et al, JHEP 10, 014 (2009), [0904.0979]. 
[32] D.-C. Dai et ai, Phys. Rev. D77, 076007 (2008), 

[0711.3012]. 
[33] D.-C. Dai et al, 0902.3577. 
[34] The ATLAS, G. Aad et al, 0901.0512. 
[35] H. Yoshino and V. S. Rychkov, Phys. Rev. D71, 104028 

(2005), [hep-th/0503171]. 
[36] U. Sperhake, V. Cardoso, F. Pretorius, E. Berti and 

J. A. Gonzalez, Phys. Rev. Lett. 101, 161101 (2008), 

[0806.1738]. 

[37] M. Shibata, H. Okawa and T. Yamamoto, Phys. Rev. 

D78, 101501 (2008), [0810.4735]. 
[38] U. Sperhake et al, Phys. Rev. Lett. 103, 131102 (2009), 

[0907.1252]. 

[39] M. Alcubierre, Introduction to 3+1 numerical relativ- 
itj/International series of monographs on physics (Oxford 
University Press, Oxford, 2008). 

[40] T. Regge and J. A. Wheeler, Phys. Rev. 108, 1063 
(1957). 

[41] F. J. ZeriUi, Phys. Rev. Lett. 24, 737 (1970). 

[42] E. Newman and R. Penrose, J. Math. Phys. 3, 566 (1962). 

[43] S. A. Teukolsky, Astrophys. J. 185, 635 (1973). 

[44] A. Coley, R. Milson, V. Pravda and A. Pravdova, Class. 

Quant. Grav. 21, L35 (2004), [gr-qc/0401008]. 
[45] R. Milson, A. Coley, V. Pravda and A. Pravdova, Int. J. 

Geom. Meth. Mod. Phys. 2, 41 (2005), [gr-qc/0401010]. 
[46] V. Pravda, A. Pravdova, A. Coley and R. Milson, Class. 

Quant. Grav. 21, 2873 (2004), [gr-qc/0401013] . 
[47] M. Ortaggio, V. Pravda and A. Pravdova, Class. Quant. 

Grav. 24, 1657 (2007), [gr-qc/0701150]. 
[48] F. R. Tangherlini, Nuovo Cim. 27, 636 (1963). 
[49] H. Kodama and A. Ishibashi, Prog. Theor. Phys. 110, 

701 (2003), [hep-th/0305147]. 
[50] U. Sperhake, Phys. Rev. D76, 104015 (2007), [gr- 

qc/0606079]. 

[51] M. Shibata and T. Nakamura, Phys. Rev. D52, 5428 
(1995). 

[52] T. W. Baumgarte and S. L. Shapiro, Phys. Rev. D59, 

024007 (1998), [gr-qc/9810065]. 
[53] H. Kodama, A. Ishibashi and O. Seto, Phys. Rev. D62, 



064022 (2000), [hep-th/0004160]. 
[54] E. Berti, M. Cavaglia and L. Gualtieri, Phys. Rev. D69, 

124011 (2004), [hep-th/0309203]. 
[55] V. Cardoso and J. P. S. Lemos, Phys. Lett. B538, 1 

(2002), [gr-qc/0202019]. 
[56] H. Yoshino, T. Shiromizu and M. Shibata, Phys. Rev. 

D72, 084020 (2005), [gr-qc/0508063] . 
[57] U. Sperhake, B. J. Kelly, P. Laguna, K. L. Smith and 

E. Schnetter, Phys. Rev. D71, 124042 (2005), [gr- 

qc/0503071]. 

[58] U. Sperhake et al, Phys. Rev. D78, 064069 (2008), 
[0710.3823]. 

[59] Cactus Computational Toolkit, 

http: //www. cactuscode . org/ 
[60] E. Schnetter, S. H. Hawley and I. Hawke, Class. Quant. 

Grav. 21, 1465 (2004), [gr-qc/0310042]. 
[61] Mesh refinement with Carpet, 

http: //www. carpet code . org/ 
[62] J. Thornburg, Phys. Rev. D54, 4899 (1996), [gr- 

qc/9508014]. 

[63] J. Thornburg, Class. Quant. Grav. 21, 743 (2004), [gr- 
qc/0306056]. 

[64] M. Ansorg, B. Bruegmann and W. Tichy, Phys. Rev. 

D70, 064011 (2004), [gr-qc/0404056] . 
[65] M. Boyle and A. H. Mroue, Phys. Rev. D80, 124045 

(2009)^ [0905.3177]. 
[66] E. Berti, V. Cardoso and A. O. Starinets, Class. Quant. 

Grav. 26, 163001 (2009), [0905.2975]. 
[67] E. Berti, V. Cardoso and C. M. WiU, Phys. Rev. D73, 

064030 (2006), [gr-qc/0512160]. 
[68] E. Berti, V. Cardoso, J. A. Gonzalez and U. Sperhake, 

Phys. Rev. D75, 124017 (2007), [gr-qc/0701086]. 
[69] P. Anninos, D. Hobill, E. Seidel, L. Smarr and W.-M. 

Suen, Phys. Rev. D52, 2044 (1995), [gr-qc/9408041]. 
[70] G. Lovelace et al, 0907.0869. 

[71] V. Cardoso, J. P. S. Lemos and S. Yoshida, JHEP 12, 

041 (2003), [hep-th/0311260]. 
[72] B. S. DeWitt and R. W. Brehme, Ann. Phys. 9, 220 

(1960). 

[73] R. H. Price, Phys. Rev. D5, 2419 (1972). 

[74] E. S. C. Ching, P. T. Leung, W. M. Suen and K. Young, 

Phys. Rev. D52, 2118 (1995), [gr-qc/9507035]. 
[75] V. Cardoso, S. Yoshida, O. J. C. Dias and J. P. S. Lemos, 

Phys. Rev. D68, 061503 (2003), [hep-th/0307122]. 
[76] S. W. Hawking, Phys. Rev. Lett. 26, 1344 (1971). 
[77] M. Davis, R. Ruffini, W. H. Press and R. H. Price, Phys. 

Rev. Lett. 27, 1466 (1971). 
[78] P. Anninos, D. Hobill, E. Seidel, L. Smarr and W.-M. 

Suen, Phys. Rev. Lett. 71, 2851 (1993), [gr-qc/9309016]. 
[79] M. E. Araujo and S. R. Oliveira, Phys. Rev. D52, 816 

(1995) . 

[80] O. M. Moreschi and S. Dain, Phys. Rev. D53, R1745 

(1996) , [gr-qc/0203071]. 

[81] D. V. Gal'tsov, G. Kofinas, P. Spirin and T. N. Tomaras, 
Phys. Lett. B683, 331 (2010), [0908.0675]. 

[82] D. V. Galtsov, G. Kofinas, P. Spirin and T. N. Tomaras, 
JHEP 05, 055 (2010), [1003.2982]. 



